Noname manuscript No. 

(will be inserted by the editor) 



Naoto Shiraishi 



Anomalous System Size Dependence of Large 
Deviation Functions for Local Empirical Measure 



Received: date / Accepted: date 



Abstract We study the large deviation function for the empirical measure (the time-averaged density) 
of diffusing particles at one fixed position. We find that the large deviation function exhibits anomalous 
system size dependence in systems with translational symmetry if and only if they satisfy the following 
conditions: (i) there exists no macroscopic flow, and (ii) their space dimension is one or two. We 
investigate this anomaly by using a contraction principle. We also analyze the relation between this 
anomaly and the so-called long-time tail behavior on the basis of phenomenological arguments. 
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I Introduction 

The large deviation theory has been intensively studied in non-equilibrium statistical physics [H-TT2]. 
In particular, since the discovery of the fluctuation theorem, which claims a symmetric property of 
the large deviation function for the time-averaged entropy production [T], large deviation functions 
for various time-averaged quantities have been investigated. These include the additivity principle 
in non-equilibrium steady states [2][3] , some relations in non-equilibrium thermodynamics 4, 5j , the 
glass transition [5], the photon emission [7j, frequency of earthquakes [8], and an operational method 
for calculating the time-averaged current [S] . As stated above, the importance of the large deviation 
functions for time-averaged quantities has been increasingly recognized. 

Here, let us consider the scenario of observing time-averaged quantities in laboratory experiments. 
In many cases, we measure such quantities at one fixed position because measurements at one position 
is much easier than those over a large area. Therefore, many studies on large deviation functions for 
such time-averaged "local quantities" arise in various contexts such as the molecular motors in living 
organisms [TO], transports in mesoscopic systems [TTJ, and stochastic processes in applied mathemat- 
ics [12]. With the above mentioned background, in this paper, we study a new aspect of the large 
deviation function for time-averaged local quantities. 

Specifically, we investigate the large deviation function for the occupation frequency at one fixed 
position, which we refer to as a local empirical measure, in systems with translational symmetry. The 
local empirical measure at a position s 6 R d is defined as 

Mas(r) ■=- I p x {t)dt, (1) 
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Table 1 Anomaly in the variance of fi in the long-time limit and the thermodynamic limit. 



dimensions 


taking the long-time limit first 


taking the thermodynamic limit first 


one 


= 0{L) 


= e (i/V7) 


two 


= 9(bxL) 


= 0(1ht/t) 


three and more 


= 0(1) (standard dependence) 


= & (1/t) (standard dependence) 



where p x (t) is a fluctuating density at x and at time t, and d is the space dimension. The large 
deviation function for H x (t) is denoted by /(/x). That is, 

ProbOu x (T) = /x)~e-^H (2) 

Here, "A ~ B" indicates lim T _ ) . 00 In A/t = lim T _ ) . 00 In B/t. The main discovery of this paper is the 
anomalous system size dependence of I(fJ>). 

Naively, we expect that I(fJ-) does not depend on the system size in the thermodynamic limit because 
I(fJ,) is a local quantity. For systems in more than three dimensions, such a naive expectation holds 
true. Surprisingly, however, the naive expectation fails for some systems in one and two dimensions. 
Concretely, in the thermodynamic limit, the variance of /i x (t), which can be calculated from I(fi), 
diverges in proportion to system size L in one dimension, and In L in two dimensions. 

This anomaly is also understood as the problem that the order of the long-time limit and the 
thermodynamic limit cannot be exchanged (see Table [1]). If we fix r and take the thermodynamic 
limit L — > oo first, Prob(/i x (T) = /n) does not satisfy the formula @, which is the definition of the 
large deviation function. Specifically, the variance of H x (t) decays as &(\/^Jt) in one dimension and 
as (9(ln t/t) in two dimensions for r — > oo, which is slower than (9(1/t) observed in standard cases. 
Here, is defined as follows: f(x) = 0(g(x)) indicates that there exist positive real numbers a and b 
such that a ■ g(x) < f(x) < b ■ g(x) for any x larger than some number. Although the existence of the 
anomalous time dependence is already known in some specific models such as random walks [13 and 
symmetric simple exclusion process (SSEP) |14j . its universal features have not been demonstrated. In 
this paper, we present general arguments for the singularity in connection with the anomalous system 
size dependence. 

This paper is organized as follows. In Sec. II, for the purpose of grasping properties of the anomaly, 
we analyze a solvable microscopic model, specifically independent random walks on a lattice with 
discrete translational symmetry. By solving this model, we suggest the essential conditions for this 
anomaly. In Sec. HI, for systems whose distribution functions follow the Fokker-Planck equation, we 
derive this anomaly by using a contraction principle. In Sec. IV, from a phenomenological viewpoint, 
we analyze the relation between this anomaly and the long-time tail behavior [15H20] . 



II Microscopic model - a solvable example 

In this section, in order to confirm the existence of the anomaly, we study a simple solvable example. 
Consider a d-dimensional cubic lattice with a periodic boundary condition. The length of the lattice in 
direction r £ {x : y, z, ■ ■ ■ } is L r . The site in the lattice is denoted by i — {i r } € N d , where < % r < L r . 
Suppose that there are N particles on the lattice. A position of the s-th particle at time t is denoted 
by x s (t). 

The particles do not interact with each other. The time evolution of each particle is described by 
a discrete-time Markov chain with a transition matrix S, where 

S j>i =Proh(x s (t + l)=j\x s (t)=i). (3) 

The migration length of a particle in one transition is bounded by a finite value independent of L r . 
Namely, by setting a r = i r — j r mod L r with —L r /2 < a r < L r /2, a necessary condition for Sjj ^ 
is that \a r \ is less than a given finite value independent of L r in each direction r. We also assume that 
the transition matrix S is periodic with a period of M r in each direction r. For example, when d = 2 
and the periods of transition matrix in x and y direction are M x and M y (and there are integers l x 
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and l y , which satisfy M x l x = L x and M y l y — L y ), Sj } i — Sj> y i' holds if i — j = i — j and i x = i' x 
mod M x , i y = i' y mod M y . 

We denote by P(x(0)) the initial probability distribution. Now, we define the local empirical mea- 
sure at the origin as 

^ TV r 

Mo(V) := - ^2^2 S o, Xs (t)> ( 4 ) 

s=l t=0 

where 5 is the Kronecker delta. First, we take the long-time limit r —5- oo. Then, we take N and 
l r to be sufficiently large under the condition that the particle density p :— NJ\\ L r , M r , and 
Sj i(0 5= V < M r ) are fixed. In this setting, we analyze the system size dependence of the variance of 
Mo(t). 

By using the fact that particles move independently, we derive the variance of Po(t) as follows. 
We denote the local empirical measure of one particle by Vq{t) : = l/ T Et=o^o^i(t)- Calculating the 
variance of vq{t) and multiplying this value by N — Y[ r L r p, we obtain the variance of /x (r). 



II (A) Systems in one dimension 



Let M be the period of the transition matrix S, and L (= Ml) be the length of the system. We define 
a transfer matrix A(h) as 

A jti (h) := e hS °.*S j>i} (5) 

where < i, j < L — 1 and ft, is a real number. We denote by [x(i)]£_ a path for the particle 1 from 
time t = to t = r. We also denote by P([x(t)]l =0 ) the probability for paths with initial probability 
distribution P(x(0)). 

We define () as an ensemble average of trajectories generated by the transition matrix S. Then, 
(e k " ( T )) satisfies 

(e hTV °^ := e,tT</o(T)p (M*)]r=o) 

[z(t)]I=o 

t=T-l 

= Yi e hTMr) n s x{t+1)Mt) p(x(o)) 

= eA(h) T P Q . (6) 

^2[x(t)] T represents the summation over all possible paths. Pq is a vector representation of P(x(0)). 
e indicates an L-dimensional vector (1, 1, 1 • • • , 1). Using the transfer matrix method, we obtain 



lim - log (eA(h) T P ) = log A max (h), 

t— >oo T 



(7) 



where A max (h) represents the maximum eigenvalue of the matrix A(h). Note that maa; (ft<), the corre- 
sponding right eigenvector to A max (h), and Po are not orthogonal. It is because Pq is a nonnegative 
vector and 4> max {h) is a positive vector, which is implied by the Perron- Frobenius theorem. It follows 
from Eqs. © and © that 



Jm, r"->„lr)"). - ( ^ ) log .l u> ,,,.(/,) 



(8) 



h=Q 



where (} c represents cumulants. We define cumulants as X^^=o h n /n\-(x n ) — exp (J^^Lo h n /n\ ■ (x") c ) . 
The first cumulant (vq(t)) c is equal to the expectation value, and the second cumulant (vq{t) 2 ) c is 
equal to variance. 

Now, our goal is to derive A^a^h). The matrix A{h) is written as A(h) — S+ (h + h? /2J X +0(h 3 ), 
where X is defined as 

/ Soft 5b, l . . . «So,Mi-l \ 
... 



X := 



\ 



(9) 
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Let v € Q L be the right eigenvectors of S. It follows from Bloch's theorem [21] that v is written as 

h 



k e*T<> 

(10) 

koe^ 6 



where 9 = 2-Kb, b = -1/2 + 1, -1/2 + 2 • • • Z/2 when Z is even, and b = -(I - l)/2, -(/ - l)/2 + 
1 • • • (Z — l)/2 when I is odd. ko ■ ■ ■ /cm-i are coefficients. Note that (ko ■ • • /t'Af-i) T is an eigenvector 
of a matrix S' G K MxM , where S' jti := EL=o 

_ Sj+ X M,i+yM and T represents a transposition. 
We can easily check that (Til)]) are eigenvectors of S and that (flT)]) takes M x Z kinds of different 
vectors. Thus, all eigenvectors of S is written as (fT0| . Let A{, i?Tl be the m-th largest eigenvalue of S 
when 9 = 2nb, and Vb, m and u^.m be the corresponding right and left eigenvectors. Vb, m and ut,^ m are 
normalized as Ub t m v b,m = 1. We denote the value of kj in i>b,m as kj(b, m). The maximum eigenvalue 
of S is Ao,o = 1- Thus, A max (h) and (p max (h) are expanded in /i as 



A max (h) =l+(h + 



h 2 



ai 



2 



0,2 



4>max( h ) = V 0,0 



h + Y 



b,tn 
(b,m)^(0,0) 



Pb,m^b,m 



b.m 



Qb,mVb,m + 



(ii) 

(12) 



where oi, <i2, P6,mj and g{, jm are coefficients. By setting n — 2 in Eq. (|8|). the variance of v§(j) is 
expressed as 

lim t (vo(t) 2 ) = 2ci2 + ai — af . 

7" — too C 

Here, a\ and 02 are written as 

i-i M-l 

ai = "0,0^0,0 = ThX] X] /c i(°'°) S 'o,A«+j, 

»=0 j=0 

_ 1 \ -> UofiXVb^m 1 \ - 



WO.O-XVft,' 



0,0 



b,m 



1 - A 



-A„ 



b.m 



(13) 

(14) 
(15) 



(fc,m)^(0,0) (b,m)/(0,0) 

iio,o^' u fc,m is in proportion to 1/ Ml. Note that ai = 0(1/1). When Z is sufficiently large, we obtain 



2a 2 c± B 



1 

Ml 



E 



A 



b.m 



1 - Ai, 



(16) 



(6,m)/(0,0) 



where B is a constant independent of Z. The dominant contribution to the right hand side in (JTSJ) is 
by terms with large Xb, m / (1 — ^b,m)- F° r large Z, the eigenvalue Xb, m takes a value near 1 only when 
m — and |6/Z| <C 1. Therefore, abbreviating A^o to Aj,, we can rewrite the right hand side in (Til)]) as 



2a 2 ^ 5 



1 

MZ 



6^0 



A 6 



(17) 
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Here, because we can reselect a larger M if necessary, without loss of generality, we assume that 
particles cannot move more than M sites in one transition. Then, we define a matrix 



S b 



( So,a 

■ 2tvI 



S . 



1 + Sm,i£ 
Si,i 



Sm.m-iz 1 1 



S\ M-l + Sm+1,M-1G 



2 lt& 

I 

2nl, 



(18) 



\S. 



M— 1,0 



'A/ 



-l,M e ' "M-1,1 + "JAf— l,Af+ie 1 



S 



M-1,M-1 



The maximum eigenvalue of 5& is X b , and let the corresponding right and left eigenvector be ip b and 
£ b . ij> b satisfies ip b :— (ko(b, 0), • • • , kM-i(b, 0)) . When \b/l\ <C 1, Sf, is expanded as 



/ 2nb 1 /2tt6 



Z + O 



(19) 



where Z is given by 



Z := 



/ 

<5l,M 



-Sm,i 




— Sm,m-i \ 

-<S f Af+l,M-l 



(20) 







\Sm-i,m 5m-i,m+i •• 

Here, Z corresponds to an operator that gives a particle current passing from Xm to 2?m-i- A& is 
expanded as 

'2tt6 x 2 
T 



2vr6 t 
Ab = 1 + i— — ■ J 



C 



(21) 



Here, C is a quadratic coefficient, which depends only on Z and So- J is defined as J := £ Zip and 
J indicates particle flow in steady states. 

Eq. (IT7|) takes qualitatively different values depending on whether J = or J ^ 0. For J = 0, 
which means that there exists no macroscopic flow, the right hand side in (|17p is 



\Ml) fel 



A h 



1 
MZ 



b#0 



/ 2 



4tt 2 & 2 • C 



D 



(22) 



for large /, where D is a constant independent of I. Here, we use the fact that the dominant contribution 
to the right hand side in (|17|) are terms with small b and other terms make a negligible contribution. 
We also use the equality X^b^i V^ 2 = ?r 2 /6. Finally, by using (TT^l) . (^ ( T ) 2 )c satisfies 



(Mo(r) 2 ) c = i>p<^(T) 2 ) o = e(L). 



(23) 



This result shows that the variance of /xo(r) diverges in proportion to L in the thermodynamic limit 
L — > co. 

Conversely, in the case J ^ 0, the right hand side in (fTT|) is 



Ml) ? 1 - A* ~ ( MZ 

• Z7T 



A 



MO 



1 



1 + iZsk . j 



6^0 



J+( 2 f) -c 



-) y 



c — J 



Ml J ^ J 2 

M0 



(24) 



for large Here, we use the fact that J^mo^ 111 ^ 6 / 0- ~ ^i>) = 0' which is implied by ImA^ = — ImA_j 
and Re Aj = Re A_j. It can be seen from Eq. (j2"4"]l that 



</,o(r) 2 > c = Lp^o(r) 2 > c -0(l), 



(25) 



where 0(1) indicates that the left hand side in (f2"5]l is independent of L. This result shows that the 
variance of Ho(t) converges in thermodynamic limit when there is macroscopic flow. 
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II (B) Systems in two and more dimensions 



For systems in more than two dimensions, we can perform calculation in a manner similar to that in 
one dimension. Therefore, in this subsection, we show an outline of the calculation. 

We consider a two-dimensional lattice. Let M x and M y be periods of transition matrix in directions 
x and y, respectively, and L x = M x l x and L y = M y l y be the length of lattice. € L x xL y represents 
a position on the lattice. Let v € Q L * xL v be the right eigenvectors of transition matrix S E R( L * xL «) . 
In a similar manner to one-dimensional case, it follows from Bloch's theorem that v is written as 



ko,M v -l 

k fle 'y y 



&M X -1,M„-1 

«Mx-i,oe ly y 

j h-- 

kM x -l,Mv-l e h 



ir-9 v 



\ko,M y -ie '* k Mx -i,M y -ie '* k 0> M y -ie ll±x(lx e h 9y ••• fcAf e -i,M„-ie r '* " x e l '* ° v J 

(26) 

where a = 2nb (b = -l x /2 + l, -l x /2 + 2- ■ ■ l x /2 when l x is even, and b = -{l x - l)/2, -(l x - l)/2 + 
1 • • • (Zjt — l)/2 when Z^ is odd) and = 27rc (c = —l v /2 + 1, —l y /2 + 2 • • • Z y /2 when Z y is even, and 
c = —(l y — l)/2, — (Z a — l)/2 + 1 • • • (Z y — l)/2 when l y is odd). fco,o • ■ • fcM«,-i,Af„-i are coefficients. 
Here, we write the vector u as a matrix only for convenience. Let Xb iC be the corresponding eigenvalue. 
Now, we define a transfer matrix A(h) in a manner similar to Eq. ([5]). and derive vl mQa; (/i), the maximum 
eigenvalue of through a perturbation expansion in powers of h. Eq. (|T3]) still holds, and Eq. (|17l) 

is modified as 

a^B^^—V E (27) 



• iaj — 1 a 
7 I- 2 ? 

kM x -i,oe 



M X Myl y l X 



X 



b.c 



(6,c)#(0,0) 



At c is expanded as 



^6,c — Aq o + i 



2?r6 2ttc\ ^ ^2tt6 + 2ttc 



(28) 



If J ^ 0, we can easily show that (^q{t) 2 ) c is independent of Z, and this implies the absence of the 
anomaly. Thus, let us consider the case when J = 0. In this case, the second term of the right hand 
side in (|2"8"|) is zero. For simplicity, assume l x = l y = I. By substituting Eq. (|2"8"|) to Eq. (|2"7|) . we obtain 




Here, we use Ao o — !• We multiply the above equation by a number of particles M x M y l 2 p, and obtain 

(Mr) 2 ) c = 0(lnZ). (30) 

This result means that the variance of /Uo( T ) shows a logarithmic divergence on I in the thermodynamic 
limit. 
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When a system is in three dimensions and there exists no macroscopic flow, we calculate 02 in a 
similar manner. The result is 

1 \ 2 x ^ \,c,d 



2a2 ~ B { M x M y M z P ) £ T 



^ ( M M M 73 ) 



b,c,a 
(6,c,d)^(0,0,0) 

2 



M x M v M z l 3 J ^ Att 2 (b 2 + c 2 + d 2 ) ■ C 

b,c,d 
(6,0,(2)^(0,0,0) 



1 \^ P 

I M x M y M z l 3 J 4n 2 r 2 ■ C 
B ^3 



(31) 



(M x M y M z ) I 3 
Multiplying Eq. by M x M y M z l 3 p, we find that 

(/i (r) 2 ) c = 0(1). (32) 
In other words, the variance of Hq(t) converges in the thermodynamic limit. 

In this section, we calculate the variance of Ho(t). We can calculate higher-order cumulants in 
a manner similar to that for variance, and they show the same anomalous system size dependence. 
Because the large deviation function and the cumulant generating function are transformed to each 
other via the Legendre transformation, the above results indicate that shows the anomalous 

system size dependence. 

Note that the "no macroscopic flow" condition (J = 0) is wider than the equilibrium condition. We 
can easily check that a model which breaks detailed balance but satisfies J = condition also shows 
the anomaly. On more realistic situation, for example, a system that is microscopically irreversible and 
macroscopically reversible J2] and a system with a non-Gaussian noise |23j exhibit no macroscopic 
flow but are out of equilibrium. Such systems break detailed balance only locally, and the anomaly 
occurs when detailed balance is broken globally. Therefore, we expect such systems also shows the 
anomaly. 



Ill Diffusive systems 

In this section, we consider diffusive systems whose probability distribution functions evolve according 
to the Fokker-Planck equation 

d -f± + V-j( X )=0 (33) 

with current densities 

j(x) := 0Dfi(x) ■ (F - VU{x)) - DVfx(x). (34) 

Here, D is the diffusion coefficient, U(x) is the potential, F is the force, and (3 is the inverse temper- 
ature. Note that random walks and SSEP systems, which are already known to show the anomaly of 
the local empirical measure (in Ref. |13U14j and the previous section), give diffusion equations similar 
to Eq. (|33[) in the hydrodynamic limit (25) . Therefore, analyses in this and the next sections can be 
applied to such systems. 

For such systems, C. Maes et al [J] found exact formulas of the large deviation functions for the 
empirical distribution function defined as fi(x,T) :— (fJ-x{T~)) xeM d- Note that ^ x {t) is a function of only 
t; in contrast, fi(x, r) is a function of both r and x. We can calculate the local empirical measure, 
/ i o( T ) 5 with the empirical distribution function. Let Idp(^(x)) be the large deviation function for the 
empirical distribution function fj,(x,r), namely 

Prob (fj,(x, t) = n{x)) - e -^WM(*0). (35) 
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Hereafter, we abbreviate /j,(x) to Ref. [4] showed that 7df(m) can be written as 

o-(^) - oy(/x) 

JdfW = 1 , (36) 

where o~(/z) is the entropy production with fi, V is a modified potential under which /i is a stationary 
solution of Eq. (|33[) . and oy(/i) is the entropy production in the modified system with V. In other 
words, when /j, is given, V is determined such that 

V • j v (x) = V • (pDfi(x) ■ (F - VV(as)) - DVfi(x)) = (37) 

is satisfied for all x. 

In this section, we calculate the local empirical measure by using a contraction principle [24] 

/(/i ) = inf / D f(mO»0) (38) 

,u(as) 

under constraint conditions /x(0) = /io and J fi(x)dx = 1. Although we can consider systems with 
any number of dimensions, in order to avoid some technical difficulties, we focus on one-dimensional 
systems in the argument below. Henceforth, for simplicity, we set U = 0. 



Ill (A) Properties of empirical distribution function 

We consider a one-dimensional periodic boundary system with length L, and take a coordinate — L/2 < 
x < L/2. We use a contraction principle to the origin x = 0. Note that the meaning of L in this section 
is different from that in the previous section. In the previous section, L is a length before taking the 
hydrodynamic limit. In contrast, in this section, L is a length after taking the hydrodynamic limit. 
However, in terms of the L dependence, the discussions in this and previous sections give the same 
results. 

The entropy production is written as 

,00 = ,3 (/ Fjdx - | / U, d x) + *M = J ^ dx , (39) 

where s(fi) denotes the Shannon entropy s(/x) := J filn^idx. Here, the first term in (|3Uf represents 
heat dissipation caused by the force, the second term in ([5^]) represents heat dissipation caused by the 
change of the potential, and the third term in (|39[) represents the change of the Shannon entropy. It 
follows from Eq. (J34J) that 

a(a) = L(3 2 DF 2 fi+ [ D^^-dx, (40) 
J M 

cvin) = L(3Fj v , (41) 
where /2 represents the uniform measure. Eq. Q34p is transformed to 

3v ■ - PD(F VV{x)) D^M, (42) 

where we set U = V. By integrating Eq. (|4"2")) around the ring, j v satisfies 

C L/2 dx 

jv — = L/3DF. (43) 

J -L/2 V 
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Set 5fj, := fj, — [A, and assume that 5/i is much smaller than /i. It is implied by Eq. (1361) that 



i / r L/2 D(Vfj,y 



dx + L 2 f3 2 DF 



'-L/2 M 

r L/2 pjy^f 

l-L/2 M 

f L ' 2 D{Vtf 



dx + L 2 f3 2 DF 



dx + (3 2 DF 2 fi 2 



l-L/2 



r L/2 D{Vnf 

'-L/2 




dx + p 2 DF 2 p 2 (I % + % + O (Sfi 3 )) ~dx 

J-L/2 \M MM / M , 



1 f L/2 (D{V5n) 2 I3 2 DF 2 5h 2 



dx. 



-L/2 



(44) 



L/2 

In the last line, we have used the relation J_ L i 2 S/j.dx = 0. 



L/2 ' 



III (B) Contraction 

By using the contraction principle in relation to Eq. (144[) . we obtain a minimizing problem 



I(Ho) = inf - 

dfj, 4 



1 f L/2 (D(V5^) 2 2 DF 2 b l i 2 



dx 



-L/2 



with constraint conditions 



dSfj, 
dx 

L/2 
-L/2 



5/i(0) = 6fx , 
d5ii 
dx 



x=-L/2 



5fj,(x)dx = 0, 



x=L/2 



(45) 

(46) 
(47) 

(48) 



where we set 5fXo := fiQ — \i. Let 6fi s (x) be a function that minimizes the right hand side in (|45p . By 
noting the constraint condition (I48|) . we use the variational method and obtain a differential equation 



V^S/jl s (x) - f3 2 F 2 5(i s (x) = const. 



(49) 



Eq. (|49j) has qualitatively different solutions depending on whether F = or F ^ 0. When f = Ono 
flow case), the solution to Eq. (14^1) is derived as 



Sfla(x) = < 



L 2 
L 2 



L 



x — 



L 

X+ 2 



for < x < #, 

2 - - 2- 



(50) 



2 



for - f < x < 0. 



Substituting Eq. ([50|) into Eq. (g5]), /(Mo) satisfies 



/(Mo) 



4/7. 



L 2 



12 



= O 



(51) 
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This result shows anomalous dependence on L. With respect to variance, Eq. (fSTj) indicates that the 
variance of fio is proportion to L, and this is consistent with results in Sec. II. 
In contrast, when F^O (flow case), the solution to Eq. (|^|) is expressed as 

_ j A + Be- Sx + Ce Sx for < x < § , 
dMs ~ \A + Be Sx + Ce- S * for - § < z < 0, ( J 

where 5 = |/3F|. The constraint conditions (gBJ), 63), and (g5]) determine A, £?, and C as 

| (^ S -l)^o 

1 + LS/2 + (LS/2 - 1) e~ LS ' 1 j 

£ = g 5 x (54) 

^l^ + ll^e-' (55) 
Substituting Eq. (JS5J) into Eq. (J3SJ), we find that 

7(/i ) = 0(1). (56) 
This result shows that there exists no anomalous dependence on L when flow exists in the system. 



IV Phenomenological arguments 



In the previous sections, it is suggested that the necessary and sufficient conditions for this anomaly 
are i) there is no macroscopic flow, and ii) the dimension of systems is one or two. In this section, we 
derive these conditions from a phenomenological viewpoint. 

We consider fluid systems that show the long-time tail behavior of density. The long-time tail is a 
phenomenon in which time correlation functions of Ap, a fluctuation of density, decay not exponentially 
but as power law. This phenomenon was first found in simulations of hard spheres as a property of 
velocity autocorrelation functions [T3]. Thereafter, it has been found that many conserved quantities 
in various fluid systems show such phenomena [16H20) . 

We assume the specific form of the long-time tail behavior as follows. Let Apo(t) be a fluctuation 
of density at the origin x = and at time t. We consider a d-dimensional cubic system with sides of 
L and suppose that there exists no flow. In this setting, we assume that the density autocorrelation 
function satisfies 



(Ap o {0)A Po (t)) { 



= /(*) 
oc t- d ' 2 

oc c~* 



< t < C, 
C<t< g(L), 

g(L) < t, 



(57) 



asymptotically in large t. Here, C denotes the time when the long-time tail behavior starts to appear, 
and C is independent of L. f(t) represents the behavior of autocorrelation functions until t = C. 
g(L) satisfies g(L) oc L 2 . This relation is true if the system shows normal diffusion. If the system size 
is infinite, (Ap o (0)Ap o (t)) oc £~ d/2 for C < t. The symbol "< ^-'V 2 " indicates that a speed of 
convergence to under t — > oo is the same as or faster than t~ x ~ d / 2 . 

Conversely, in the case that the system size is infinite and there is macroscopic flow, we assume the 
relation 



(Ap o (0)Ap o (t)) < 



= /(*) 
oc t- d ' 2 

< t -l-d/2 



< t < c, 
C <t< h(F), 
h{F) < t. 



(58) 



Here, F is the force that causes the macroscopic flow. h(F) is a function of F, and h(F) satisfies 
h(F) oc F~ 2 . The symbol "< t~ 1 ~ d / 2 " indicates that a speed of convergence to under t — > oo is the 
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same as or faster than i~ 1_d / 2 . Note that the formula of f(t) and the value of C arc different from Eq. 
(fSTjl . Systems considered in Sec. HI show such long-time tail behavior ((57)) and (see Appendix). 

A mathematical structure of the anomaly of the local empirical measure is very similar to that of 
the long-time tail. Owing to this, in this section, with the assumption of the long-time tail (|57[) and 
(|58p. we derive the anomalous dependence of (/xo(r) 2 ) c on both r and L. 



IV (A) Derivation of anomalous time dependence 



Suppose that there exists no flow in a system. First, we take the thermodynamic limit L — > oo, and 
then, we take a large r. By using Eq. (|57|) . we calculate the dependence of (pq(t) 2 ) c on r as 



<Mo(r) 2 ) c = ( 



1 



JO 



db'dt"Apo{t')Ap Q {t")) 



-i • 2 / T di"'(r - t'")(/\/5 O (0)Zipo(t'")> 

T J 



2 It 



dt"'f{t"')+ / tft"'(T - t'") ■ at"'- d / 2 



'± -2(tK + a'T 2 - d ' 2 ) ford ^2, 
4 • 2 (rif + o't In r) for d = 2, 



(59) 



where K = J dtf(t). a and a' are constants. In the third line, we have used an approximation 
r-C~r. 

It follows from Eq. (JSHJ) that the dependence of (po( T ) 2 )c on T is 

/- 1 



(Mo(t) 2 ) c oc < 



Inr 

T 



for d = 1, 



for d = 2. 



for d > 3. 



(60) 



This result represents the anomalous dependence on r in one and two dimensions. 



IV (B) Derivation of anomalous system size dependence 



Suppose that there exists no flow in a system. First, we take the long-time limit r — > oo, and then, we 
take a large L. By using Eq. (|57p . we calculate the dependence of (po( T ) 2 )c on L as 

lim (M T ) 2 }< 



lim • 2 / df"(r - H<^Po(0)^po(r')> 



lim 



2 T 



s(i) 



dt"'f{t"')+ / tft'"(T - *"') • at'"- d / 2 



' Urn i • 2 frKT + a'r f g(L) x ~ d l 2 - C 1 ^ 2 )) for d ^ 2, 



lim — 

r— )-oo 7"* 



• 2 rif + a'r In 



for d = 2, 



(61) 



where a and a' are constants. It follows from Eq. (j6"Tj) that the dependence of (ho(t) 2 ) c on L is 



lim (a*o(t) 2 ) c oc < 



L for d = 1, 
lni ford = 2, 
const for d > 3, 



(62) 



for large L. This result represents the anomalous dependence on L in one and two dimensions. 
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IV (C) No anomaly with macroscopic flow 

Suppose that macroscopic flow exists in a system. First, we take the thermodynamic limit L — >■ 00, 
and then, we take large r. From Eq. (|5T|) . the dependence of (po( T ) 2 )c on r is calculated as 

(Mr) 2 ) c = (\ f T f dt'dt"Ap (t')Ap (t")) 



J 



= ^ • 2 f dt"'(r - f ")<2Vpo(0)^po (*"')> 
T •/ 

1 / fC MF) ,r \ 

~ _L . 2 T / dt"7(i'") + r / dt"'bt"'- d ' 2 + / dt'"(r - *'") • ai'"- 1 -^ 2 

T \ Jo Jc Jh(F) J 

~ i • 2 (rif + a'r 1 "^ 2 ) , (63) 



where a, 6, and a' are constants. As implied by d > 1, Eq. ((63]) shows standard dependence on t; 



(^o(t) 2 ) c a -. (64) 

T 



This result shows that the anomaly vanishes when flow exists in a system. 



V Concluding remarks 

The essential conditions for the anomaly in systems with translational symmetry are i) there exists no 
macroscopic flow, and ii) the dimension of the system is one or two. When these conditions arc satisfied, 
a power law behavior of statistical properties is observed. Concretely, in Sec. IE, we show that /x s (:r), 
the function minimizing I(p(x)), decays from the origin x — as a power law in the space direction. 
In Sec. IV, we show that the time correlation function of Ap (t) has a power law tail structure in the 
time direction, and the anomaly is derived directly from this structure. Hence, this anomaly can be 
seen not only for the local empirical measure but also for other time-averaged conserved quantities 
that show the long-time tail behavior, such as the velocity [T5HT5] and the stress tensor |20j . 

In particular, for systems in which total momentum is conserved, we expect that variance for 
time-averaged pressure shows anomalous time and system size dependence even in three dimensions. 
Concretely, when we take the long-time limit first, the variance behaves as O(L), and when we take 
the thermodynamic limit, the variance behaves as @(l/yr) for any dimensions. The reason why the 
anomaly remains even in three and more dimensions is as follows; for a system in d-dimensions, a 
side wall is in d — 1-dimensions. The difference between the dimensions of systems and side walls is 
always one, and thus, with projections, the variance for time-averaged pressure seems to show the same 
anomaly as that for the one-dimensional case. This is a future problem. 

In conclusion, we have two comments. First, in Sec. HI we show that the anomaly is caused by 
the contraction process. Here, it is known that some nonequilibrium steady systems show anomalous 
fluctuations called long-range correlations [26]. When this anomaly occurs, a fluctuation decays not 
exponentially but as power law in the space direction. Recent studies show that the long-range corre- 
lation can be understood as the anomaly caused by a contraction process for the time direction, and 
the derivation also uses a modified potential 3 . This is very similar to the discussion in Sec. IE, and 
therefore the long-range correlation and the anomaly we discuss in this paper seem to have a common 
origin. This is also a future problem. 

Second, in recent years, new universal relations in nonequilibrium thermodynamics have been in- 
vestigated [27, 28]. With regard to this development, Ref. |S] has proposed a novel type of inequality 
in nonequilibrium thermodynamics, in the sense that it is written with the large deviation function 
for the empirical distribution function. It would be remarkable if the anomaly reported in this paper 
is related to nonequilibrium thermodynamics. 

Acknowledgements The author thanks S.-I. Sasa for providing useful advices and fruitful discussions. The 
author also thanks M. Otsuki and T. Nemoto for providing helpful comments. 
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Appendix: Proof of the long-time tail behavior 

In this appendix, we prove the long-time tail behavior (|5T[) and (|58|l for diffusive systems whose 
distribution functions evolve according to the Fokker-Planck equation ((33|) with [7 = 0. Suppose that 
the system is in d dimensions under periodic boundary conditions with length L. We also assume that 
the initial condition is Ap x (0) — for x ^= 0. 

First, for F = 0, we derive Eq. (|57p . The evolution of Ap is described by 

d -^ = D^(A P ). (A.1) 
From the Fourier transforms of Eq. (|A.1[) , it is shown that 

zww = E e ^ |fe|2t > ( A - 2 ) 

k 

where k takes k = 2n/L ■ j for all j € N d except j = 0. When 4Dir 2 t/L 2 < 1, we can transform Eq. 
(|A.2I) to an integral form as 

Ap (t) = j e- D ^ 2t dk = t~i J e- D \ k '\ 2 dk' cc t~ d / 2 . (A.3) 

Here, we set y/ik = k! . In contrast, when AD-K 2 t/L 2 ^> 1, the right hand side in (|A.2[) is written as 

oo 

Y^e- D W 2t ~ £ ^ D|fe|2 ' + E E e-^l^d-e-^l^occ-*. (A.4) 

fe |fc|=27r/L ra=2|fc|=27iVn/£ 

Here, we use 

OC OO OO , 2 ^ 

E E e ^' fe|2i ^ E( 2d )" e_4 ^ = E e(-^ +ln2d )" = o (e- D m 2t ) (A.5) 

"=2 \k\=2TT^/n/L n=2 n=2 

Eqs. (|A.3|) and (|A.4[) satisfy the long-time tail behavior (|57l) . 

Next, for inhnite size systems with F ^ 0, we derive Eq. (1551) . The evolution of Z\p is described by 

^p- = -pDF-V{Ap) + DV 2 {Ap). (A.6) 
From the Fourier transforms of Eq. (|A.6[) , it is shown that 

Ap (t) = J e(- D ^ 2 -^ DF - k >dk = J e- D ^ 2t cos(I3DF ■ kt)dk. (A.7) 
When (3D \ F\ t <C y/Dt, we perform calculation in a manner similar to Eq. (|A.3|) as 

Aoo(t) = y e-^l 2 ' cos((3DF ■ kt)dk ~ / e'^'dk = t~i j e~ D \ k '\ 2 dk' oc f~ d / 2 . (A.8) 



In contrast, when /3D |f | i ^> v-Di, we use an approximation for large A 



k 2 i rk 2 ,2 f 

/(fc) cos(2^4fc)dfc ~ — ^ ^dfc, (A.9) 



which follows from 

r' k '+i i d 2 f 

(A.10) 
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We divide the vector k into fcj_ and fe|| (fe = k± + k\\), where k± is perpendicular to F, and k\\ is 
parallel to F. By using the approximation (|A.9|) . we obtain 

J e- DWH cos(/3DF ■ kt)dk = J j er D{\k ± \ 2 +k\)t cos(/3jD | F | kllt ) dk±dkll 



= C-t-^-- t / [-D + D 2 k\t) e- Dfc i*dfc|| 



<t^ d/2 . (A.ll) 
Here, C is a constant. Eqs. (IA.8I) and (lA.llj) satisfy the long-time tail behavior |58|) . 
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